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Abstract 

We present an improved version of our earlier work on summing 
the planar graphs in (f)'^ field theory. The present treatment is also 
based on our world sheet formalism and the mean field approximation, 
but it makes use of no further approximations. We derive a set of 
equations between the expectation values of the world sheet fields, and 
we investigate them in certain limits. We show that the equations can 
give rise to (metastable) string forming solutions. 



*This work was supported in part by tlic Director, OfRce of Science, Office of High 
Energy and Nuclear Physics, of the U.S. Department of Energy under Contract DE-AC03- 
76SF00098, in part by the National Science Foundation Grant 22386- f3067-44-X-PHHXM, 
and in part by the Department of Energy under Grant No. DE-FG0f-97ER-4f 029 



1. Introduction 

One of the most challenging problems in the study of non-Abelian gauge 
theories is to discover the dual string description [1]. After many years of 
relative stagnation, the discovery of the AdS/CFT correspondence [2,3] has 
led to important progress in the resolution of this problem. This approach 
relies heavily on ideas from string theory and supersymmetry, and it can deal 
effectively with only a fairly restricted class of field theories. 

Recently, we initiated a program [4], one of the goals of which is to estab- 
lish the duality between field and string theories. The idea, which goes back 
to Nielsen, Olesen, Sakita and Virasoro [5], and which was given a systematic 
foundation via 't Hooft's 1/Nc expansion [6], is to construct a world sheet 
description of the planar graphs of a general field theory. As shown in [6], 
planar graphs are selected by introducing an internal color degree of freedom 
and taking Nc ^ oo limit, and the use of a mixed coordinate- momentum 
space hght cone coordinates leads to the world sheet description. Although 
originally only a (j)^ theory was considered for simplicity [4, 6], this approach 
is sufficiently flexible to be able to handle more reahstic cases, including 
nonsupersymmetric [7] and supersymmetric [8] gauge theories. 

The original world sheet model suggested by 't Hooft was non-local. In 
reference [4], it was shown how to reformulate this model as a local world 
sheet field theory by introducing additional non-dynamical fields. The world 
sheet structure of a given graph consists of the bulk and a bunch of "solid" 
lines, which describe the multiple internal boundaries representing the loops 
of a multi-loop planar diagram. The dynamics is all in the boundaries, and 
the only function of the fields that live in the bulk is to instantaneously 
transmit the interaction from one boundary to another. Therefore, the local 
world sheet theory can be thought of as a topological theory. 

The representation of planar graphs by a topological field theory is per- 
haps elegant, but it has no new physical content. The elimination of the bulk 
fields would lead back to the original non-local structure of 't Hooft. One 
gets something new and interesting if the topological bulk fields are some- 
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how promoted into genuine dynamical fields. One way for this to happen is 
through the condensation of the solid fines (the boundaries) . In a graph of a 
given order, the density of the solid lines, by which we mean the percentage 
of the area on the world sheet occupied by solid lines, is zero. This is because 
solid fines are fines; they have zero thickness. However, as the order of the 
graph asymptotes to infinity, solid lines become more and more numerous 
and dense, and ultimately, one can envisage a limit in which they acquire a 
finite density on the world sheet. This is what we mean by the condensa- 
tion of the boundaries. In this limit, the distinction between the bulk and 
the boundary disappears, the world sheet acquires a uniform texture, and 
it becomes possible to have a string formulation of the sum of the Feynman 
graphs in terms of dynamical fields on the world sheet. We will call this 
mechanism string formation through condensation of boundaries. Whether 
this really happens depends on the dynamics: One clearly needs an attractive 
interaction and also an effectively strong coupling, leading to the domination 
of higher order graphs. 

The worldsheet formalism of [4] might also provide a natural setting for 
understanding confinement in real QCD, where strong coupling is only a fea- 
ture of infrared dynamics. Then the physical world sheet representing the 
confining fiux tube would not have the uniform texture mentioned above 
but would include important fluctuations in which regions of the world sheet 
would not have solid line condensation and would retain the topological char- 
acter of perturbation theory. These fluctuations would describe the point-like 
structures implied by asymptotic freedom. 

In an earlier paper [9] , using 0^ theory as a prototype toy model, we inves- 
tigated the possibility of boundary condensation leading to string formation. 
Although this is an unphysical theory, it provides a simple setting for devel- 
oping the tools needed to attack more interesting and also more complicated 
theories. The main tool used in [9] was the mean field approximation, or the 
self consistent field approximation. This approximation has been widely used 
in both field theory and many body theory, in most cases leading to at least 
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qualitatively reasonable answers. For the problem at hand, the easiest way 
to implement it systematically is to consider the limit of large D, where D is 
the number of transverse dimensions in the light cone picture. We should add 
that although the large D limit provides a convenient bookkeeping device, it 
does not truly capture the physics of the problem in most cases. The critical 
factor is the number of degrees of freedom of the system in question: The 
mean field approximation is usually successful when applied to a system with 
a large number of degrees of freedom. Of course, one way to to have a large 
number of degrees of freedom is to have a large number of space dimensions, 
but this is rather artificial in most cases. A more physical situation is to have 
a large number of particles in a many body system, and this is the case in 
most of the standard applications. In our case, as explained above, we are 
interested in graphs with a large number of boundaries, and we can roughly 
identify the number of boundaries with the degrees of freedom. With this 
identification, the mean field method should be a good tool for investigating 
problems involving graphs with a large number of boundaries. There is then 
no need to appeal to the large D limit, except, as a convenient bookkeeping 
device. 

In this article we aim to improve the treatment of [9] in several respects. 
First, we applied the mean field approximation to a system of Ising spins 
introduced as a representation of the sum over all arrangements of solid lines. 
Since the world sheet on which these spins lived remained 2 dimensional even 
as — > oo, it was not clear why that should justify replacing each spin by 
a mean field. Instead the large D hmit is actually that of an 0{D) vector 
model. As is well-known this limit justifies treating the 0{D) invariants, the 
scalar products of the vector fields, classically. In fact the process of replacing 
the scalars with classical fields eventually leads to a mean field treatment of 
the Ising spins as well, but, unlike in [9], there are no ambiguities in how to set 
up the action for the classical mean fields. We incorporate this clarification 
in the current article. But the most substantial improvement we make here 
is in the treatment of the Dirichlet boundary conditions. In [9] we introduced 
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Gaussian representations of delta functions of the type 



2n 



) 



D/2 



5(q) = lim 




to enforce the boundary conditions. Keeping (3 finite amounts to impos- 
ing an infrared cutoff on target space. With a sufficiently exact calculation 
there is nothing wrong with this. But after making an approximation to the 
cutoff theory, the delicacy of removing the infrared cutoff might well invali- 
date the approximation. In the current article we shall impose the boundary 
conditions exactly, thereby eliminating this difficulty. The price will be the 
need to introduce a more complicated system of mean fields, some of which 
are nonlocal on the world sheet. The presence of such nonlocal mean fields 
might even make it easier to incorporate the non-homogeneous worldsheet 
textures required for a realistic worldsheet description of QCD. We are then 
able to derive a set of equations for these fields, which relies on no other 
approximation except for the mean field approximation. We consider these 
equations (eqs.(46) and (47)) as the main result of this article. Although 
these equations are non-linear and somewhat complicated, it should be pos- 
sible to attack them by various approximation schemes and numerical meth- 
ods. However, in this article, we will be content with studying them in some 
limiting cases, and we leave a full investigation of this problem to future 
research. 

The paper is organized as follows. In the next section, we present a brief 
review of the world sheet formalism developed in [4] for summing planar 
graphs. We also argue that, in the leading order of the mean field approx- 
imation, a good deal of simphfication takes place: The ghost fields can be 
dropped, and the world sheet dynamics is described by a free action (eq.(7)), 
plus constraints on the boundaries given by eq.(8). In section 3, the boundary 
constraints are rc-cxprcsscd in a form suitable for inclusion in the effective 
action and also for taking the large D limit. This is done by introducing 
auxiliary scalar fields, which arc equal to the scalar products of the vector 
fields in D dimensions. The important point is that these scalar fields become 
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classical in the large D limit. Section 4 is devoted to the construction of the 
effective action as a function of the scalar fields. This is done by explicitly 
carrying out the functional integration over the vector fields. The classical 
equations for the scalar fields resulting from this effective action completely 
determines the large D dynamics. In section 5, we investigate these equa- 
tions in various limits, and we argue that one solution, although not stable, 
can lead to condensation of the boundaries and string formation. This is a 
moderately encouraging but preliminary result, although no definite conclu- 
sions can be reached until we have the full solution. Section 6 summarizes 
our conclusions. In the appendix, we show that the contribution of the ghost 
fields is unimportant in the leading order of the large D limit. 

2. A Brief Review 

We will be working with a massless (jr" matrix field theory in the large 
N limit, which amounts to a summation of the planar graphs. As 't Hooft 
has shown [6], Feynman rules are especially simple if a particular mixture 
of coordinate and momentum light cone variables are used. We will use 
the following notation: A Minkowski vector will be written as (f f ~, v), 
where = {v^rLv^) / \/2^ and the boldface letters label the components along 
the transverse directions. The Lorentz invariant product of two vectors v and 
w is given hy v -w — y - w — v'^w~ — v~w'^. The evolution parameter (time) 
is x'^, and the Hamiltonian conjugate to this time is p~. A massless on- 
shell particle thus has the energy p~ = p^/2p+. The number of transverse 
dimensions D is arbitrary to start with, and eventually, we will consider 
the large D limit. As explained in [4], this limit is a convenient method of 
organizing the mean field approximation. Of course, the eventual case of 
interest is D — 2. 

Let us now briefly review the Feynman rules derived in [6]. A specific 
graph is represented by a set of parallel solid lines drawn on the world sheet. 
A propagator corresponds to a strip bounded by two solid hues (Fig.(l)). If 
p is the momentum carried by the propagator, the strip has width p+, and 
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Figure 1: A Propagator 

length T — x'^ . We associate two transverse momenta qi and q2 with the 
two sohd hnes forming the boundary. The transverse momentum p of the 
propagator is their difference: 



p = qi - q2, 

and using a Euchdean world sheet metric, the propagator is given by 



2p 



+ 



exp (-T (qi - qs) V V) , 



(1) 



Now let us consider a more complicated graph, with interaction vertices, 
pictured in Fig. (2). Interaction takes place at the points where a solid line 
ends, and a factor of g is associated with each such vertex, where g is the 
coupling constant. After putting together the propagators and coupling con- 
stants associated with the graph, one has to integrate over the positions of 
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Figure 2: A Typical Graph 



the sohd hnes and the position of the interaction vertices, as well as the 
momenta carried by solid lines. We note that momentum conservation is 
automatic in this formalism. 

These Feynman rules look quite non-local; however, it was shown in ref- 
erence [4] that they can be reproduced by a local field theory defined on the 
world sheet. Later, this result was extended from 0^ to more complicated 
and more reahstic theories [7, 8] . The world sheet field theory in question 
can be formulated either on a latticized or continuum world sheet, and here 
we choose the continuum version. We will eventually need to specify some 
cutoffs, but we defer that until the need arises. Using the coordinates a in 
the direction and r itself in the — t direction, the bosonic fields q((T, r) 
and the fermionic fields 6((T, t) and c((T, t) , called ghosts, are introduced. In 
contrast to q, which has D components, h and c each have D/2 components. 
(Here, we are assuming that D is even). The free part of the action on the 
world sheet is given by 

-So = £ dr da (b' ■ c' - ^q'^) , (2) 



where the prime denotes derivative with respect to a. Here p"*" is the + com- 
ponent of the total momentum flowing through the graph. For convenience, 
we have also restricted the r integration; eventually, we will let — > — oo 
and Tf +00. This action is to be supplemented by the following boundary 
conditions: On solid lines, q((7, r) is constrained to be independent of r, or, 
equivalently, the Dirichlet boundary condition 

is imposed, where the dot denotes derivative with respect to r. Also, if the 
transverse momentum carried by the whole graph is p, the constraint 

r^daq' = p (3) 

Jo 

has to be imposed. For simplicity, we will set p = in what follows, which 
means periodic boundary conditions on q: 

q((7 = 0)=q((7=p+). (4) 

The corresponding boundary conditions on the ghosts are simple: 

6 = c = (5) 

on the sohd lines. 

In addition to 5*0, it was shown in [4] that the full action contains interac- 
tion terms in the form of the insertion of ghost vertices. These insertions are 
needed to correctly reproduce the factor l/2p'^ in front of the exponential 
in the expression (1) for the propagator. If, for example, one integrates the 
matter and ghost fields of Sq with the boundary conditions appropriate for 
the propagator, without any vertex ghost insertions, the result is 

^(r)exp(-rpV2p+), (6) 

instead of (1). Local ghost vertex insertions can be arranged to correct this 
defect. However, we will now argue that this is the correct leading term in 
the large D expansion. To see this, let us rewrite the propagator in the form 

e{T) exp (-T pV2p+ - ln(2p+)) . 
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The first term in the exponent really corresponds to D terms since it is the 
square of a vector with D components. Consequently, in the large D limit, its 
contribution is of the order oi D. In contrast, the contribution of the second 
term is D independent. Therefore, in the leading order in D, which is all we 
are going to consider in this paper, the factor of l/2p+ can be dropped. This 
means that, in the large D limit, we are entitled to work with Sq given by 
eq.(2) and ignore the ghost insertion terms. 

We are going to take advantage of one more simplification. Although 
the fermionic (ghost) part of 5*0 does contribute in the large D limit, it will 
be shown in the appendix that its contribution harmlessly shifts some of 
the fields we had already introduced. So, bearing in mind that the fields 
might need to be interpreted differently, we are going to drop the ghost 
contributions. 

To summarize, the starting point of the present work will be the action 



on solid lines and the periodicity condition given by eq.(4). We would like 
to stress that the terms we have neglected do contribute in the non-leading 
orders of the large D expansion. 

3. The Boundary Conditions 

Our goal is to express the sum over all the planar graphs in the form of an 
effective action. We start with eqs.(7) and (8), which generate the simphfied 
Feynman rules of the large D limit. We would like to rewrite the conditions 
(8) in a form that can be identified as a contribution to the action, and 
cast them in a form suitable for taking the large D limit. We note that the 
boundary conditions decouple in the a direction, and so it is convenient first 
to discretize the a coordinate into small segments of length Act = ei, and 





plus the Dirichlet boundary conditions 



q=0 



(8) 
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to impose the conditions for each discrete value of a separately. (In [9] this 
cutoff ei was called m the minimal unit of p"*" in the discretized formulation.) 
The parameter ei will play the role of an ultraviolet cutoff on the worldsheet. 
It should be pointed out that we are making an important change in the 
way we classify graphs. The customary classification is according to the 
total number of solid Uncs, which is the same as counting the powers of the 
coupling constant. This is analogous to a Fock space description. Instead, 
we are now going to focus on the distribution of solid lines separately for each 
value of a. This is more analogous to the occupation number description. 
If one is interested in the condensation of the solid lines, the occupation 
number description is clearly superior, since counting the number of lines in 
a condensate is not very useful. 

Next, for a fixed value of o", we are going to carry out the summation 
over all possible partitionings of solid lines and over all values of momenta 
fiowing through them. This will result in an integral equation. In particular, 
we note that the Dirichlet boundary conditions (8) will be taken into account 
exactly. In this respect, the present treatment is superior to the one given 
in [9]; in that work, only a weaker version of these boundary conditions was 
imposed. 

In Fig. (3), we have drawn lines located at the discrete values of a and 
extending in the r direction. These lines consist of alternating solid and 
dotted segments: As before, the solid lines are located where the boundary 
conditions (8) hold, whereas the dotted lines run through what used to be 
blank space, where there are no boundary conditions. Now consider a typical 
line located at some o", starting at and ending at a variable point r. For 
simplicity, we assume that this line contains an equal number of solid and 
dotted line segments, and we call this number m. One has then to sum 
over all possible partitions of this line between alternating solid and dotted 
segments for a given m, and then sum over all m. Let Fjn{a,Ti,T) be the 
factor that takes care of (8) on the solid segments. This function satisfies the 
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Figure 3: Solid and Dotted Lines 



following recursion relation: 

JTi Jti 

In going from to F^+i, we have added a dotted segment extending from ti 
to T2, and a solid segment extending from T2 to r. Since there is no boundary 
condition on the dotted segment, we associate with it a factor of one. With 
the solid segment, we associate the factor K{a,T2,T), which is introduced 
to take care of the corresponding boundary condition. The total F is then 
given by 

00 

F{a,n,T)= J2Fm{<y,n,T). (10) 

m=0 

and it satisfies the integral equation 

F{a, Ti, t) = Fo{a, Ti, t) + dn / dT2 K{a, T2, T)F{a, n, ri). (11) 

So far, we have defined an F that implements (8) for a single line located 
at a. To find the F that implements these boundary conditions for the whole 
world sheet, we have to multiply the F's associated with each discrete value 
of a and let r — > t/: 

Fin,Tf)=l[F{ar,,T„Tf). (12) 
n 
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Now let us return to the integral equation (11). We will set 

F^{a,T,,T) = 5{T,-T). (13) 

The kernel which has to implement (8) on the solid line segment extending 
from r2 to T, is given by 

K{(T, T2, t)^ j dkj D\ exp (i £ dr' \{a, r') ■ (q((7, r') - k)) . (14) 

The integration over 1((T, r') sets q(o", r') equal to a constant r independent 
vector k in the interval {t2, r), and then k is integrated over. It is clear that 
this is equivalent to the condition that q is independent of r' in this interval. 

It will become clear as we proceed that an ultraviolet cutoff in the variable 
k is needed to avoid divergences. We therefore modify the expression for K 
by introducing a suitable cutoff of the form of exp(— e/c^) and then do the 
integration over k: 

K{a, Ti, T2) ^ I dk 1 L>1 exp (^-e^ + i dr'l{a, r') ■ (q((7, r') - k)) 

~Te £ /r ■ ^"0 • ^^^^ 

We note that in addition to the worldsheet ultraviolet cutoff ei introduced 
earlier, we now have a second ultraviolet cutoff e, acting in the target space. 
Later, we will find that a third cutoff may be necessary, although it is likely 
that the third cutoff is expressible in terms of the first two. 

It is worth noting that the nonlocal term in eq.(15) is entirely due to 
the integration over k, the Dirichlet value of q on the internal boundary 
which K adds to the worldsheet. But in the worldsheet formalism applied 
to gauge theories with extended supersymmetry [8], the worldsheet variables 
q((T, r) were 2 + d dimensional where = 2 for A/" = 2 supersymmetry and 
(i = 6 for A/" = 4 supersymmetry. The process of dimensional reduction, 
which leads to the 4 dimensional theory with extended supersymmetry, was 
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achieved in [8] by imposing the boundary condition = on all boundaries, 
external and internal, whenever I — 3, 4, ... d. This means that the same 
components of the vector k are put to zero and not integrated. In other 
words the nonlocality is not introduced for the components of 1 in the extra 
dimensions. This worldsheet locality for the dimensions associated with the 
field theoretic internal degrees of freedom could be associated with the generic 
appearance of Sn factors of the space-time manifold on the strong coupling 
side of the AdS/CFT correspondence^ 

Wc are dealing here with a vector model, and the standard first step 
in taking the large N (in this case, the large D) limit, is first to express 
everything in terms of the scalar products of the vectors. From eq.(15), we 
see that this goal has been achieved: The kernel K is expressed in terms of 
the scalar products q ■ 1 and 1-1. In the next step, the scalar products are 
replaced by their vacuum expectation values, which are proportional to D, 
and which we denote by D 0i and D 02 respectively: 

1(<7, r') • q(<7, r') ^ {l(a, r') ■ q(a, r')) = D <t>,{a, r'), 
l(a,T')-l(a,r") ^ r') ■ l(c7, r")) = 02(a, r', r"). (16) 

We would like to make it clear that 0i and 02 are fixed classical background 
fields, not to be integrated over. They are the saddle points that domi- 
nate the large D limit, to be determined by minimizing the effective action. 
The general rule is that vector valued fields are quantum mechanical, and 
they are to be integrated over. In contrast, scalar valued fields are classical 
background fields. 

The kernel K can now be written solely in terms of 0i and 02- 

{n/ef/^ exp (^iD dr'M'^, r') - £ J^' dr' dr''^'^, r', r")) ■ 

^ ' " " (17) 

'^We thank J. Maldacena for an illuminating conversation that inspired this suggestion. 



13 



Substituting this expression back into eq.(ll), we have to solve an integral 
equation to find F{Ti,Tf). We will reconsider this problem in the next section, 
after the simplifications following from translation invariance are taken into 
account. Here we note that once F is known, the total action S, which 
incorporates the boundary condition (8) on solid lines, is given by 

exp{S)=exp{S,)F{n,Tf), (18) 

where, 

S, = £\a I ciT(-^(q'((7,T))^ + zAi(a,T)(l((j,T) •q((j,T) -L'0i((7,t)) 
- ij^dajdrj dr'A2(a, r, r') r) • r') - D 02 (a, r, r')) • (19) 

The first term on the right is 5*0 (eq.(7)). In the rest of the terms in 5"!, the 
Lagrange multipliers Ai and A2 enforce the definition of 0i^2 given by (16). 

4. The Effective Action 

In this section, we will address the problem of the computation of the 
effective action. This will involve 

a) solving the integral equation (11) for F and 

b) doing the functional integrations over the vector-valued fields q and I, 

keeping all the scalar valued fields fixed. 

We will be able to make only partial progress in solving eq.(ll). In con- 
trast, the functional integrations over q and I can be carried out in closed 
form. Once the effective action is derived, the classical fields 0i,2 and Ai,2 
are determined by the equations 

The calculations outlined above are greatly simplified by making use of 
translation invariances in both the a and r directions on the world sheet. 
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This symmetry means that various background fields are either constants or 
functions of a single variable: 



where Ai and 0i are independent of a and r, and A2 and 02 depend only on 
the difference r — r' . Substituting this in eq.(20), we see that S\ will be non- 
local in the r coordinate. This makes the problem of finding the minimum of 
the effective action more difficult; instead of an algebraic equation, we have 
to solve functional equations in one variable. 

Let us now go back to the integral equation (11), taking advantage of the 
simplification resulting from translation invariance. First, we note that the 
functions K and F simplify: 



Now, eq.(ll) can be formally solved using Fourier transforms. We define 



The lower limits of integration start at zero, since both F and K are defined 
to be zero for negative values of the argument. The solution to eq.(ll) is 
then given by 



Ai(a,r)^Ai, 0i((T,r)^0i, 
A2 (cr, T, r') ^ A2 (r - r') , 02 ((J, r, r') ^ 02 (r - r') , 



(21) 



K{(y, Tl, T2) K{t2 - Ti), F{a, Ti, T2) F{t2 - Ti). 



(22) 




(23) 




(24) 



where K can be written in the form 




(25) 



and where we have defined 




(26) 
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The contribution of F to the effective action follows from eq.(18): 

S = S, + Sf, SF = hi{F{Tf-Ti)). (27) 

Now, let us examine Sp in more detail. The time interval Tf — Ti, which 
corresponds to r in eq.(23), will eventually tend to infinity. It is a standard 
result in Fourier transform that this limit will be dominated by the lowest 
lying singularity in the variable cu of the integrand in eq.(23), which could 
be the tip of a cut or a pole. Since, being conjugate to time, cu is energy, 
a cut corresponds to the continuum and the pole to a bound state. We 
are particularly interested in bound states, since they can lower the energy 
of the system and lead to a non-trivial solution to our variational problem. 
Therefore, we will assume that the lowest lying singularity is a pole*, and we 
will be looking for the zeroes of the denominator in eq.(24): 

u - 27iig^K{cu) = 0. (28) 

If cuo is a solution to this equation, the corresponding energy is given by 

Eo = iuo, (29) 

and the corresponding contribution to the action is 

SF^--p+{Tf-Ti)uJo, (30) 
^1 

where ei is the cutoff resulting from the discretization of the interval along 

the a direction. (See eq.(12)). At this point, we should note that bound 

states cannot arise in perturbation theory, so our treatment from now on is 

definitely non-perturbative. There is also the question of the reality of Eq. 

In a stable theory, the spectrum should be real and bounded from below; 

however, since we are dealing with 0^, an intrinsically unstable theory, we 

■tin reference [9], the (uiergy spectrum consisted of only discrete states due to the 
infrared cutoff introduced in the Gaussian representation of the delta functions imposing 
Dirichlet boundary conditions. 
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may end up with complex energies. We will investigate this question in 
section 5 in the context of a simple model. 

For the time being, we cannot proceed any further with the bound state 
problem without knowing 02, so we will leave its solution in the implicit form 
defined by eqs.(24), (25) and (28). Let us now turn to the next problem, 
that of carrying out the functional integrals over the quantum fields q and 1. 
Replacing Ai and A2 in cq.(19) by their translation invariant forms (eq.(21)), 
this equation can be rewritten as 

Si = s[ + s';, 

S[ = iD da J cir Ai0i + J dn J cir2 A2(ti - T2)02(ti - T2)^ , 

S'l = £ da j dr (^-^q!^ + iXil ■ 

- dajdrij (iT2 A2(t2 - Ti) 1((7, Ti) • 1((7, T2). (31) 

It turns out to be more convenient to work in the momentum space. We 
define 

q((T,T) = / dpoY^ex.p{ipia + ipor) q{po,pi), 
^ PI 

A2(t) = j dpo ex.p{ipoT)X2{po): 

M^) = j dpo exp{-ipoT)(f)2{po), (32) 

where 

Pi — 27rm/p'*", m & Z, 

because of the periodic boundary conditions in the a direction. 

Before attempting to do the functional integrals over q and 1, we will 
examine the classical equations of motion. They will enable us to make a 
comparison with the standard string action. Written in momentum space. 
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the equations of motion are 



iXil — p\q = 0, 
Aig-47rA2(po)^~ = 0. 



(33) 



The solution to these equations is 




(34) 



The usual string action corresponds to taking 




(35) 



where a is the slope parameter. One can see this as follows: In coordinate 
space, eq.(35) is equivalent to the action 



The Hamiltonian corresponding to this action is the light cone energy p , 
and the quantized values of p" are given by the usual string result 

mr 

P = —r, 
ap^ 

where n is a positive integer. We are really interested in the squares of the 
masses of the excitations, and these are given by 



since we have set (see eq.(4)) p = 0. Prom this result, we see that the 
slope parameter is given by a = {2Tq)~^ = ira', where a' is the slope of open 
string Regge trajectories. A string that is linearly confining at large distances 
corresponds to a A2 that has the dependence given by the eq.(35) for small 
Po, so we are going to investigate whether such a behavior is compatible with 
our dynamical scheme. 
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After this diversion, let us now carry out the integrations over q and I in 
S'{. The result is the determinant, which is the product of the eigenvalues k, 
that satisfy the equations 

iXil — p\q = Kq, 
iXiq — 4:niX2{po)l' = i^^i, (36) 

or 

+ (pI + 47riA2(po)) + + 4:niplX2{po) = 0. (37) 

The eigenvalues depend on both po = P, which is continuous, and on 

27rm 

Pi = —r 
p+ 

with m an integer. After the integration over q and I, we end up with the 
factor 

(det)-°/^ = nn U-i^pM + A?)"" 

m p \ P ) 

The corresponding Tr In is given by 

(38) 

In writing the right hand side of this equation, wc have factored A2 from the 
first term on the right, and we incorporated it into the second term in the 
form of ln(A2). The coefficient of this term involves a sum over m and appears 
to be divergent. However, remembering that the a variable was discretized 
into segments of length ei, this sum is easily seen to be equal to p^ /^i- 
We now collect various contributions to the effective action: 

S,^S[ + SF--^Tr\n. (39) 

The terms that appear in this formula are given in eqs.(31), (30) and (38). 
It is now easy to write down various classical equations of motion that follow 
from this action. The equation obtained by varying with respect to A2 is 

2^^Dp\r, - r^^,{p) - ^§2^ = 0, (40) 

2 d{X2{p)) 
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or 

= .-^i— coth ( ^s'/\p)] 1 , (41) 

"^'^"^^ (47r)3Ai(p)sV2(p) [2 Snh,X2{p) ^ ' 

where, for convenience, we have defined 



s{p) 



47rA2(p) 



Eq.(41) is the first fundamental equation of our dynamical scheme; it 
estabhshes a relation between 02 and A2. A second relation follows from 
varying the action with respect to 02 : 

In this case, it turns out to be more convenient to stay in position space and 
use L instead of 02 as the independent function. The connection between 
them is 

202(t) = L"(t), 

where we have assumed that 02 (p) is an even function of p. This assumption, 
although not crucial, simplifies matters. It is consistent with all our equations 
and with the approximate solutions given in Section 5. 

Taking into account the contribution of S[ (eq.(31)) and of Sp (eq.(30)), 
we have 

We have now to compute the functional derivative of uoq with respect to L 
from eq.(28). Differentiating this equation gives 

(43) 

and from eq.(25), we find 

6iK{uJo)) D (, DL{t)\ 

7(LMr = ~^W exp(^..or + ^i^0.r-^J. (44) 
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Eqs(42),(43) and (44) put together give us the second relation between 
02, or equivalently L, and A2: 



^„^^^ ^ igoexp[iu;oT + zD(^,T- ^ ^^^^ 

2eei (1 + dr t exp (iujoT + iDcpiT - ^^^^^ ^ ' 

where the dimensionless couphng constant gQ is defined by 



TT 



\ -0/2 



To the above equations, one has to add the equations resulting from 
varying the action with respect to Ai and 0i: 

Ml = 47r y dp (a2(p)02(p) + , 

fOO 

Ai = 2ie / dTTX'^ir). (46) 
Jo 

We now have a complete set of equations needed to solve for the classical 
fields 01,2 and Ai,2- We would like to stress that, apart from the large D 
limit, which is the basis of the mean field method, so far everything is exact. 
Unfortunately, to make progress, we need to know ujq, which is only given 
implicitly through eq.(28). Ultimately, it should not be too difficult to find 
approximate solutions by numerical methods. We will, however, leave this 
problem to future research. Instead, in the next section, we will try to extract 
as much information as we can from these equations by analytic methods. 
Such information should be useful for the eventual numerical work. 

5. Consequences of the Dynamical Equations 

In the last section, we have derived a complete set of equations for the 
dynamical variables Ai_2 and 0i,2- This section will be devoted to the inves- 
tigation of the consequences of these equations. First, for the convenience 
of the reader, we collect the the set of equations to be investigated; namely, 
eqs.(28),(41) and (45): 

roc / D \ 

= uq — ig^ J dr exjp iiujQT + iDcpiT — —L{t)j , 
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Mp) = COth i , 

■ 9 /■ , ■ n DL(t)\ 

Kir) = y " ^ y ^' ^ (47) 

2ee,{l +g'oI^drrexp{iu;or + iD(|>^T-^)) 



L(r)= / ^(2-e^^--e-^^-)02(p), (48) 



where 

dp 

which follows from eq.(26). 

We will now try to find an approximate solution to eqs.(47), and the 
approximation will be based on a pole dominance model. This is a crude 
model, which is at best valid only for large r, or equivalently, for small p. We 
remind the reader that this region is of interest in probing string formation 
(see the section following eq.(35)), so any information gleaned is of value. 
Our starting point is the ansatz 

h ^ Pp' (49) 

as p — > 0, where /3 is a constant. The basic idea is to cycle this limiting 
behavior through the dynamical equations to check its consistency. We first 
need to restate it in the position (r) space. Consider the r ^ 00 limit of 
L(t) from eq.(48). In the terms involving the exponentials, this limit can 
be deduced from the well known Fourier relation: Large r is dominated by 
small p. Since the integrand is non-singular at p = by virtue of the ansatz 
on 02, the exponential factors oscillate to zero for r — > 00, with the result 
that 

L(t) ^ Lo, (50) 
where Lq is a constant. This constant is given by 



J^Mp)- (51) 
J p 



If L{t) goes to a constant as r — >■ 00, the integral for K{uj) (eq.(25)) 
develops a pole in the variable uj. We will now make the assumption that this 
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pole is a good approximation for K for small values of cu (pole dominance). 
We therefore have 

K(u)^—(-) exp -] -. 52 

ul + DcPiooo + f = 0, (53) 



Eq.(28) now becomes 



where we have defined 

f = gle.p[-^). (54) 
The two solutions to this equation are 

--lD<p,±l{4f-D^<l>lf\ (55) 



We will investigate both solutions. 

Next, we evaluate A2 by replacing L(r) by Lq on the right hand side of 
eq.(45). As explained above, this should be a good approximation in the 
large r, small p regime. Transforming into the momentum space, we have 

T / ^ ^ f i^O + D<Pl? 1 

'^^^ 47reei ((c^o + Dcj^^Y - g^) (a;o + Dcj)^ - p) p^ 

We note that the 1/p^ dependence of A2 for small p is exactly what is needed 
for string formation (see eq.(35)). We can now substitute this result in cq.(41) 
to determine the small p behavior of 02, to see whether it is consistent with 
our initial ansatz (49). In the limit p ^ 0, we indeed find that 

— > const -p^, (57) 

establishing consistency. 
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It is also important to establish the reality properties of various variables 
within the context of our approximate solution. We will initially assume 
that 01 is pure imaginary and is real and positive. It then follows from 
eq.(56) that, at least for small p, A2 is pure imaginary. In the same small p 
limit, eq.(57) tells us that 02 is real, and therefore, from eq.(51), so is Lq. 
We are now able to verify one of our initial assumptions: A real Lq means 
that is real and positive (eq.(54)). It then follows from eq.(55) that a;^ 
are pure imaginary, and the corresponding energies E'^ (eq.(29)) are real, as 
they should be. Next, we will investigate Ai. It can be calculated from the 
second eq.(46), where A2 is given by eq.(45), with L(r) replaced by Lq. The 
result is 



which shows that Ai is real. Finally, we go back to the first equation (46): 
Since Ai and 02 are real, and A2 is pure imaginary, it follows that 0i is pure 
imaginary, verifying the remaining initial assumption. 

Prom the above analysis, we have seen that A2(p) goes hke for small p, 
which is the necessary condition for string formation. However, in addition, 
the square of the slope parameter a (eq.(35)) must be real and positive. 
Using eq.(56), this parameter is given by 



and remembering that Ai is real and 0i is pure imaginary, is positive 
for the lower sign. This means that in cq.(55) for loq^ we must also choose 
the lower sign. This is then the string forming solution, and from the same 
equation, we see that the corresponding energy is positive. However, the 
other solution, corresponding to the upper sign, has lower energy, since E^ is 
negative, so that it will actually dominate the worldsheet path integral, and 
not describe a string. The state that looks like string is at best mcta-stable 
and more likely completely unstable. This is, of course, not unexpected, 
since 0^ is an inherently unstable theory. We expect a more realistic theory. 



Ai = 



ei - (a;o + £>0i) 



2 ' 



(58) 




(59) 
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such as a non-Abelian gauge theory, to retain the string forming solution but 
to be free of the instabihty. As explained in the introduction, our goal in 
studying 0^ is to develop the tools needed to attack more realistic but also 
more complicated theories in a simpler setting. We find it encouraging that, 
despite the instability, we find a signal for string formation. 

The results discussed above were based on the pole dominance approxima- 
tion. As we have argued earlier, this is probably a reasonable approximation 
for small p. In fact, the dependence of A2 for small p, which leads to 
string formation, appears to be generic. This is easily seen from eq.(45): In 
momentum space, the double derivative on A2 turns into a factor of p^, and 
unless the right hand side of this equation vanishes accidentally at p = 0, 
A2 must have a 1/p'^ dependence. On the other hand, we cannot expect the 
pole dominance model to be valid for large p. In fact, we show below that 
the large p regime is quite different from the small p regime. 

We will try to determine the large p behavior of the fields through self 
consistency. We start with the ansatz 

Mp) b\p\, 
for p — > 00, where 6 is a constant, or, equivalently 

Mr) ^ 

as T — > 0. It then follows from 

2<^2(r) = L"(r) 

that 

L(r) ^ 461n(r), 

in the same limit. Eq.(45) then gives 

A2(r) ^ const • r^^-*-^/^) 

again, as r ^ 0. The motivation for the original ansatz was in fact to obtain 
a power behavior for A2, which translates into 

A2(p) const • \p\''^/'-^ 
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for large p in momentum space. Now, from eq.(41), we find that 

Mp) ^ const • A(p)-3/2 ^ const • \p\^/^-^f>D/2e^ 

and comparing this with the initial ansatz, we see that it is consistent if 

bD _7 

Admittedly, this is a crude analysis, but if we accept it, we have asymptotic 
limit 

^2{p) — ^ const ■ 02(p) — const ■ \p\, (60) 

for p ^ oo. This is quite different from the small p behavior 

^2{p) — const ■ <p2{p) — const -p^. (61) 

We end this section with a short discussion of the need for a third cutoff. 
The asymptotic behavior derived above shows that the integral for L (eq.(48), 
as well as the integral for Ai0i (eq.(46)), are divergent at large p. The 
integrals have to be regulated, say, by a factor of the form 

exp(-e2/), 

where 62 is a third cutoff, in addition to e and ei. This new cutoff is probably 
not independent of e. The argument goes as follows: Instead of an exponen- 
tial cutoff, let us latticize the r direction, with a lattice spacing At. Prom 
(6), we see that one lattice step At corresponds to a factor of 

exp (^-At pV2p^) • 

Comparing this with the cutoff exp(— e/c^) (eq.(15)), we can identify 

At fti ep+. (62) 

On the other hand, there is a rough relation between At and €2'- 

62 ^ (At)2, 
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and therefore, it follows that 

£2 ~ (ep+)'- (63) 

In view of this connection between e and 62, it is tempting to conjecture 
that one of them is redundant. In fact, in eqs.(46) and (47), e and ei do not 
regulate any divergent integrals; they merely appear as parameters. One can 
easily eliminate the singular dependence on these parameters by the following 
scahng of the fields: 

02 e02, A2 — , 
eei 

Ai — , 01^01, ujo^ujo. (64) 
After this scaling, as far as the equations are concerned, one can let 

e — > 0, ei — > 0, 

while keeping the dimensionless coupling constant go fixed. Let us see how 
some physical quantities behave under this scaling. As indicated above, cuq 
and therefore Eq are independent of e and ei, whereas the slope parameter 
scales as e/ei (eq.(35)). Keeping this ratio finite as each cutoff parameter goes 
to zero would avoid any singular behavior. So it appears that at least in the 
leading order of the mean field approximation, both e and ei are redundant. 
In contrast, as pointed out earlier, the cutoff 62 cannot be dispensed with; 
it is really needed to regulate integrals over p in equations (46) and (48). 
Clearly, the problems involving cutoff dependence and renormalization must 
await a better understanding of the solution to eqs.(47). 

6. Conclusions 

The main result of this article is the derivation of a set of equations for 
the sum of all planar graphs the (f)^ field theory, subject to the mean field or 
large D approximation. Apart from this single approximation, the treatment 
is exact. This is in contrast to our earlier paper [9], where the same problem 
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was considered, but also further and somewhat questionable approximations 
were made. The present treatment is therefore a definite improvement over 
the one presented in [9] . 

The equations derived in this article are somewhat involved but they ap- 
pear treatable by numerical methods. As a preliminary step, we have studied 
them in the large and small momentum regimes. The latter is especially im- 
portant since it is relevant to string formation. We find that a string forming 
solution is consistent with our equations, although this solution is unstable. 
Given that we are dealing with an unstable theory to start with, we still find 
this result encouraging. 

Much still remains to be done. A thorough numerical analysis of the 
equations derived here should be carried out. Since we do not necessarily 
expect to find interesting results for the (j)^ model, the methods developed in 
this article should be applied to more physical theories, whose world sheet 
formulations are already available [7, 8] . It is also important to go to next to 
leading order in the large D expansion, both as a check on the leading order, 
and also to investigate questions such as Lorentz invariance. 
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Appendix. Contribution of the Ghost Sector 

In this appendix, we will discuss the changes resulting from taking into 
account the ghost fields b and c. First of all, an additional factor has to be 
added to K given by (14): 

K xKg, 

where, 

Kg{a, T2, t) = J J Dcexp (^i J dT'(b{a, t') ■ b{a, t') + c{a, t') ■ c(a, r'))^ 

(65) 

Here, b and c are lagrange multipliers that enforce the boundary condi- 
tions (5) on solid lines. Moreover, there is a ghost contribution to the ac- 
tion(eq.(28)): 

S^S + Sg, 

where, 

^9 = £ j dr {b' -c' + Kiib-b- Drii) + K2{c-c- Dr/a)) • (66) 

In this expression, the first term comes from eq.(2), and in the rest of the 
terms, the lagrange multipliers Ki^2 set the expectation value oib-b and c • c 
equal to 771^2 ■ As before, this is in preparation for taking the large D limit: 
Vector valued fields b and c are traded for scalar valued fields ki^2 and 771,2- 

Next, as in eq.(22), we use translation invariance to replace the fields Ki^2 
and 771,2 by constants independent of a and r. As a result, we can set 

Kg^e^p{tD{7]i + 7]2){r-T2)), (67) 

and eq.(18) gets replaced by 

K ^ (7r/e)^/^ exp (^iD{r2 - n){4>i + ?7i + V2) - £^(t2 - Ti)) . (68) 

We see that as far as K and therefore F is concerned, all that the ghosts 
have done is to shift 0i by 771 + 772. As a result, eqs.(25) through (28) still 
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continue to be valid, with the proviso that 0i be replaced by the shifted 0i. 
These shifts do not change anything done in this paper, since we have never 
made any use of the actual value of 0i, and that is the reason why we have 
suppressed the ghost sector in the main body of the paper. We expect them 
to become important in the next non-leading order in D. 

Although we do not need it, for the sake of completeness, we will carry 
out the functional integrations over 6, c, h and c in Sg (eq.(62)). The resulting 
Tr In is very simple: 

C 

Trhig = p+{Tf - Ti) — KlKl (69) 

where C is a cutoff dependent constant. The reason for this simple answer is 
that the term b' ■ d in (62), the only term that "propagates" the ghosts, does 
not contribute to the determinant. This is an example of the non-dynamical, 
or equivalently, non- propagating ghosts discussed in [9]. In the absence of 
propagation, the determinant, or the Tr In has purely algebraic dependence 
on the k's. In any case, we see that the ghost contribution to the action is 
relatively trivial. 

References 

[1] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Nucl. Phys. B 636 
(2002) 99 [arXiv:hep-th/0204051]. 

[2] J. M. Maldacena, Adv. Theor. Math. Phys. 2 (1998) 231-252, hep- 
th/9711200. 

[3] O.Aharony, S.S. Gubser, J. Maldacena, H.Ooguri and Y.Oz, Phys. 
Rept.323 (2000) 183 [arXiv: hep-th/9905111]. 

[4] K. Bardakci and C. B. Thorn, Nucl. Phys. B 626 (2002) 287 [arXiv:hep- 
th/0110301]. 



30 



[5] H. B. Nielsen and P. Olesen, Phys. Lett. 32B (1970) 203; B. Sakita and 
M. A. Virasoro, Phys. Rev. Lett. 24 (1970) 1146. 

[6] G. 't Hooft, Nud. Phys. B72 (1974) 461. 

[7] C. B. Thorn, Nucl. Phys. B 637 (2002) 272 [arXiv:hep-th/0203167]. 

[8] S. Gudmundsson, C. B. Thorn, and T. A. Tran, arXiv:hep-th/0209102. 

[9] K. Bardakci and C. B. Thorn, arXiv:hep-th/0206205. 



31 



